b=zeros(1,6);
aE=zeros(6,6);
aI=zeros(6,6);
gama=0.25;

b(1) = 0.15791629516167136;
b(2) = 0.;
b(3) = 0.18675894052400077;
b(4) = 0.6805652953093346;
b(5) = -0.27524053099500667;
b(6)=gama;

aE(3,1) =0.221776;
aE(3,2) = 0.110224;
aE(4,1)=-0.04884659515311857;
aE(4,2)=-0.17772065232640102;
aE(4,3)=0.8465672474795197;
aE(5,1)=-0.15541685842491548;
aE(5,2)=-0.3567050098221991;
aE(5,3)=1.0587258798684427;
aE(5,4)=0.30339598837867193;
aE(6,1)=0.2014243506726763;
aE(6,2)=0.008742057842904185;
aE(6,3)=0.15993995707168115;
aE(6,4)=0.4038290605220775;
aE(6,5)=0.22606457389066084;
aE(2,1)=2*gama;

aI(3,1)=0.137776;
aI(3,2)= -0.055776;
aI(4,1)=0.14463686602698217;
aI(4,2)=-0.22393190761334475;
aI(4,3)=0.4492950415863626;
aI(5,1)=0.09825878328356477;
aI(5,2)=-0.5915442428196704;
aI(5,3)=0.8101210538282996;
aI(5,4)=0.283164405707806;
aI(2,1)=gama;
for j=2:6
   aI(j,j)=gama; 
end
for j=1:5
   aI(6,j)=b(j); 
end